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The High Resolution Fly's Eye HiRes-I detector has now been in operation in monocular mode for over six 
years. During that time span, HiRes-I has accumulated a larger exposure to Ultra-High Energy Cosmic Rays 
(UHECRs) above 10^® eV than any other experiment built to date. This presents an unprecedented opportunity 
to search for anisotropy in the arrival directions of UHECRs. We present results of a search for dipole distributions 
oriented towards major astrophysical landmarks and a search for small-scale clustering. We conclude that the 
HiRes-I data set is, in fact, consistent with an isotropic source model. 



1. Introduction 

The observation of Ultra-High Energy Cos- 
mic Rays (UHECRs) has now spanned over forty 
years. Over that period, many source models 
have been proposed to explain the origin of these 
remarkable events. In the past five years, theoret- 
ical models have been suggested that would po- 
tentially produce dipole distributions oriented to- 
wards M87 m or Centaurus A |2I3| . In addition, 
the Akeno Giant Air Shower Array (AGASA) has 
reported findings suggesting a 4% dipole-like en- 
hancement oriented towards the Galactic Center 
present in its events with energies around 10^^ eV 
This result seemed to be corroborated by 
findings published by the Fly's Eye experiment 
in 1999 that suggested the possibility of an en- 
hancement in the galactic plane also at energies 
around 10^^ eV ;5,, and also by a re-analysis of 
data from the SUGAR array that was published 
in 2001 that showed an enhancement in the 
general vicinity of the Galactic Center. 

However, both AGASA and Fly's Eye are sub- 
ject to a limiting factor; they are both located 
too far north in latitude to directly observe the 
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Galactic Center itself. The re-analysis of SUGAR 
data actually demonstrated an excess that was 
offset from the Galactic Center by 7.5° and was 
more consistent with a point source than a global 
dipole effect 0- While the current High Resolu- 
tion Fly's Eye (HiRes) experiment is subject to a 
similar limitation in sky coverage as the AGASA 
and Fly's Eye experiments, we will show that, by 
properly estimating the HiRes aperture and an- 
gular resolution, we can effectively exclude these 
dipole source models to a certain degree of sen- 
sitivity. However, we are not able to completely 
exclude the findings of AGASA or the theoretical 
predictions mentioned above. 

Additionally, over the past decade, the search 
for sources of Ultra-High Energy Cosmic Rays 
(UHECRs) has also begun to focus upon small 
scale anisotropy in event arrival directions. This 
refers to statistically significant excesses occur- 
ring at the scale of < 2.5°. The interest in 
this sort of anisotropy has largely been fueled by 
the observations of the Akeno Giant Air Shower 
Array (AGASA). In 1999 [T] and again in 2001 
jHI , the AGASA collaboration reported observing 
what eventually became seven clusters (six "dou- 
blets" and one "triplet" ) with estimated energies 
above ^ 3.8 x 10^^ eV. Several attempts that have 
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been made to ascertain the significance of tliese 
clusters returned cliance probabilities of 4 x 10~^ 
[5] to 0.08 [ini. 

By contrast, the monocular (and stereo) anal- 
yses that have been presented by the High Res- 
olution Fly's Eye (HiRes) demonstrate that the 
level of autocorrelation observed in our sample 
is completely consistent with that expected from 
background coincidences |11I12I13| . Any analy- 
sis of HiRes monocular data needs to take into 
account that the angular resolution in monocular 
mode is highly asymmetric. 

It is very difficult to compare the results of the 
HiRes monocular and AGASA analyses. They are 
very different in the way that they measure auto- 
correlation. Differences in the published energy 
spectra of the two experiments suggest an energy 
scale difference of 30% |14I15| . Additionally, the 
two experiments observe UHECRs in very differ- 
ent ways. The HiRes experiment has an energy- 
dependent aperture and an exposure with a sea- 
sonal variability These differences make it 
very difficult get an intuitive grasp of what HiRes 
should see if the AGASA claim of autocorrelation 
is justified. In order to develop this sort of intu- 
ition, we apply the same analysis to both AGASA 
and HiRes data. 

Our methods for detecting the presence of a 
dipole source model and small-scale clustering are 
based upon comparisons between the real data 
and a large quantity of events generated by our 
Monte Carlo simulation program. The simulated 
data possess the same aperture and exposure 
as the actual HiRes-I monocular data set. We 
show how the asymmetric angular resolution of 
a monocular air fluorescence detector can be ac- 
commodated in this method. 

2. The HiRes-I Monocular Data 

The data set that we consider consists of events 
that were included in the HiRes-I monocular 
spectrum measurement |14ll6j . This set con- 
tains the events observed between May 1997 
and February 2003. There were 1526 events 
with reconstructed energies above 10^*'^ eV and 
52 events with reconstructed energies above 
10^^'^ eV observed during this time period. The 
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Figure 1. The geometry of reconstruction for a 
monocular air fluorescence detector 



data set represents a cumulative exposure of ~ 
3000 km^-sr-yr at 5 x 10^^ eV. This data was sub- 
ject to a number of quality cuts that are detailed 
in the above-mentioned papers |14I16| . We pre- 
viously verifled that this data set was consistent 
with Monte Carlo predictions in many ways in- 
cluding impact parameter {Rp) distributions jl4| 
and zenith angle distributions ^Jj . For this study, 
we presumed an average atmospheric clarity jl8| . 

In order to perform anisotropy analysis on this 
subset of data, we must first parameterize the 
HiRes-I monocular angular resolution. For a 
monocular air fluorescence detector, angular res- 
olution consists of two components, the plane of 
reconstruction, that is the plane in which the 
shower is observed, and the angle if) within the 
plane of reconstruction (see figure ^1. We can 
determine the plane of reconstruction very accu- 
rately. However, the value of tp is more difficult 
to determine accurately because it is dependent 
on the precise results of the profile-constrained fit 

The HiRes-I angular resolution is therefore de- 
scribed by an elliptical, two-dimensional Gaus- 
sian distribution with the two Gaussian parame- 
ters, a^p and (ipianc, being defined by the two an- 
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Figure 2. The arrival directions of the HiRes-I monocular data with reconstructed energies above 10^^'^ eV 
events and their la angular resolution 



gular resolutions. For events with reconstructed 
energies above 10^^'^ eV, — [4.9,6.1]° and 
CTpianc[0.4, 1.5]°. In figureEl the arrival directions 
of the HiRes-I events are plotted in equatorial co- 
ordinates along with their la error ellipses. 

In order to understand the systematic uncer- 
tainty in the angular resolution estimates, we con- 
sider a comparison of estimated arrival directions 
that successfully reconstructed in both HiRes-I 
monocular mode and HiRes stereo mode. We con- 
sider all mono/stereo candidate events with esti- 
mated energies above 10^*'^ eV. In stereo mode, 
the shower detector planes of the two detectors 
are intersected, thus the geometry is much more 
precisely known and the total angular resolution 
is of order 0.6°, a number that is largely corre- 
lated to CTpiano and thus is negligible when added 
in quadrature to the larger term, a^. This al- 
lows us to perform a comparison of the angular 
resolution estimated through simulations to the 
observed angular resolution values of actual data. 
In figure |21 we show the distribution of angular 
errors for real and simulated data. The uncer- 
tainty in the slope of the ratio (figure Ot) leads 
to an 7.5% uncertainty in the angular resolution. 



3. Dipole Measurement Results 

The method for measuring the anisotropy am- 
plitude, a, for potential dipole sources is dis- 
cussed at length in Abbasi et al. [T7|. 

3.1. Analysis 

In summary, we measure the value of the 
anisotropy amplitude by the following method: 

1. We calculated the value of <cos6'> for the 
dipole function of the real data sample. 

2. We created a total of 20,000 simulated data 
samples, 1000 each for 0.1 increments of a 
from -1.0 to 1.0, each with the same number 
of events as the actual data. 

3. We constructed curves corresponding to the 
mean and standard deviation of <cos 6> of 
the dipole function for each value of a. 

4. We determined the preferred value of a and 
the 90% confidence interval of a for each 
dipole source model by referring to the in- 
tersections of the 90% confidence interval 
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Figure 3. Arrival direction error comparison be- 
tween real data (mono vs. stereo) and simulated 
data for events with estimated energies above 
j^q18.5 rpj^g solid line histogram corresponds 
to the arrival direction error distribution of the 
monocular reconstructed Monte Carlo simulated 
data. The crosses correspond to the arrival direc- 
tions error distribution observed for actual data 
by comparing the arrival directions estimated by 
the monocular and stereo reconstructions. The 
solid line in the ratio component corresponds to 
the &t y = ax + b where a = 0.000 ± 0.011 and 
b = 0.98 ±0.11. 



curves with the actual value of <cos 6> for 
the dipole function of the real data. 

The results are shown in tabled 

3.2. Conclusion 

We are able to place upper limits on the value 
of \a\ for each of our three proposed dipole 
source models. However, these limits are not 
small enough to exclude the theoretical predic- 
tions |ll2l3j . Also, they do not exclude the find- 
ings of the AGASA collaboration in terms of 
the intensity of the dipole effect that they ob- 
served or in terms of the energy considered be- 



Source 


a 


Galactic 


0.005 ±0.055 


Centaurus A 


-0.005 ±0.065 


M87 


-0.010 ±0.045 



Table 1 

Estimation of a via the value of <cos6'> for the 
dipole function. 



cause the events in the dipole effect observed 
by the AGASA detector possessed energies be- 
low 10^*-^ eV 1^. Since it appears that angular 
resolution has little impact on the measurement 
of a and we do not appear to be systematically 
limited, we conclude that the driving factor in 
making a better determination of a will simply 
be larger event samples. HiRes-I mono will con- 
tinue to have the largest cumulative aperture of 
any single detector for the next three to five years, 
thus it will continue to serve as an ever more pow- 
erful tool for constraining dipole source models. 

4. Small-Scale Clustering Results 

The method for measuring the small scale clus- 
tering is discussed at length in Abbasi et al. ^U] . 

4.1. Analysis 

We measure the degree of small-scale clustering 
by means of an autocorrelation function. It is 
calculated as follows: 

1. For each event, an arrival direction is sam- 
pled on a probabilistic basis from the error 
space defined by the angular resolution of 
the event. 

2. The opening angle is measured between the 
arrival directions of a pair of events. 

3. The cosine of the opening angle is then his- 
togrammed. 

4. The preceding steps are repeated until all 
possible pairs of the events are considered. 

5. The preceding steps are repeated until the 
error space, in the arrival direction of each 
event, is thoroughly sampled. 
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Figure 4. The autocorrelation for the HiRes-I 
events above 10^^-^ eV. 



6. The histogram is normahzed and the result- 
ing curve is the autocorrelation function. 

A well-behaved measure of the autocorrelation 
of a specific set of data is the value of <cos 9> for 
6 < 10°. This value is also a measure of the sharp- 
ness of the autocorrelation peak at cos6' — 1. 
However, this method of quantification does not 
depend on bin width and it does produce Gaus- 
sian distributions when it is applied to large num- 
bers of sets with similar degrees of autocorrela- 
tion. An additional advantage to this method is 
that by considering the continuous autocorrela- 
tion function over a specified interval, both the 
peak at the smallest values of 9 and the corre- 
sponding statistical deficit in the autocorrelation 
function at slightly higher values of 9 are taken 
into account. Thus we simultaneously measure 
both the positive and negative aspects of the au- 
tocorrelation signal. 

Using the description of the HiRcs-I monocu- 
lar angular resolution above, we then calculate 
the autocorrelation function via the method de- 
scribed above. In figure ^ we show the result 
of this calculation. For this sample, we obtain 
<cos6'>[o=ao=]= 0.99234. 

We also calculate the autocorrelation function 
for the published AGASA events We show 
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Figure 5. The autocorrelation for the published 
AGASA events |7|. 



the result in figure [S] For this sample, we obtain 
<cos6'>[oo 10°]= 0.99352. 

To study the relative sensitivity of AGASA and 
HiRes-I, we measure the value of <cos^>[qo io°] 
for multiple simulated sets with a variable num- 
ber of doublets inserted. We then construct an 
interpolation of the mean value and standard de- 
viation of < COS0 >[o°,io°] from a given num- 
ber of observed doublets for each experiment. 
This will allow us to state the number of dou- 
blets required for each experiment in order for 
the 90% confidence limit of < cos >[oo^io°] to 
be above the background value of 0.99250. In 
general, for a HiRes-I-like data set, the 90% con- 
fidence lower limit corresponds to the mean ex- 
pected background signal with the inclusion of 
6.25 doublets. For an AGASA-hke, the 90% con- 
fidence lower limit corresponds to the mean ex- 
pected background signal with the inclusion of 5.5 
doublets. This demonstrates that while AGASA 
has a slightly better ability to perceive autocor- 
relation, the sensitivity of the two experiments 
is comparable. However, the observed HiRes-I 
signal corresponds to the 90% confidence upper 
limit with the inclusion of only 3.5 doublets be- 
yond random background coincidence. 
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4.2. Conclusion 

We conclude that the HiRes-I monocular de- 
tector sees no evidence of clustering in its highest 
energy events. Furthermore, the HiRes-I monoc- 
ular data has an intrinsic sensitivity to global au- 
tocorrelation such that we can claim at the 90% 
confidence level that there can be no more than 
3.5 doublets above that which would be expected 
by background coincidence in the HiRes-I monoc- 
ular data set above 10^^'^ eV. From this result, we 
can then derive, with a 90% confidence level, that 
no more than 13% of the observed HiRes-I events 
could be sharing common arrival directions. This 
data set is comparable to the sensitivity of the 
reported AGASA data set if one assumes that 
there is indeed a 30% energy scale difference be- 
tween the two experiments. It should be empha- 
sized that this conclusion pertains only to point 
sources of the sort claimed by the AGASA col- 
laboration. Furthermore, because a measure of 
autocorrelation makes no assumption of the un- 
derlying astrophysical mechanism that results in 
clustering phenomena, we cannot claim that the 
HiRes monocular analysis and the AGASA analy- 
sis are inconsistent beyond a specified confidence 
level. 
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